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' Abstract 

! The recent high statistics NMC data on the Tin to Carbon structure function ratio 

seems to indicate, for the first time, a significant Q 2 dependence, especially at small 
Q_i' values of Bjorken x, x < 0.05, and Q 2 > 1 GeV 2 . Pure leading twist, perturbative 

QCD-based predictions, which are consistent with the free nucleon data, yield a fairly 
fiat ratio with little or no Q 2 dependence. In view of this seeming contradiction, we 
re-examine the applicability of such a perturbative model to nuclear structure functions 
in such a kinematical regime. We find that the model is consistent with all data, within 
experimental errors, without any need for introducing additional higher twist contri- 
butions. The model correctly reproduces the Q 2 dependence of the Carbon structure 
\Q . function as well. We also critically examine the Q 2 dependence of the corresponding 

| spin dependent structure functions. 

G\ 
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"Sh! 1 Introduction 

Q.. 

The structure function, F^(x,Q 2 ), for various nuclei of mass number, A, such as 4 He, 6 Li, 
12 C, 40 Ca, and 208 Pb, has been precisely measured in several lepton-nucleon deep inelas- 
tic scattering (DIS) experiments [jl], 0, over a large range of the Bjorken variable, x, and 
momentum transfer, Q 2 . In all cases, the ratio of the structure function for the bound nu- 
cleon in a nucleus, A, to that for the free nucleon, D (taken to be an "average" nucleon in 
Deuteron), R A = Ff/Ff, was found to be less than unity at small and large values of x, 
with antishadowing (R > 1) at intermediate x values around x = 0.1. No significant Q 2 
dependence of these ratios has been observed so far jl], |], within the error bars of these 
measurements. That is, the nuclear (bound nucleon) and free nucleon structure functions 
seem to exhibit the same Q 2 dependence. 

A recent measurement || by the NMC presents, for the first time, a high-statistics study 
of the Q 2 behaviour of the Tin to Carbon structure function ratio, F 2 Sn /F 2 c ; consequently, 
the statistical errors are much smaller than in the earlier measurements. The data appear 
to show significant Q 2 dependences, unlike what has been observed before, for other nuclei. 
Several theoretical models of nuclear structure functions are available 0, EL H], that can 
account for such a Q 2 dependence, especially at small values of x and Q 2 . Many of these 
use a hadronic approach to the problem, for example, vector meson dominance (VMD), 
which is valid down to nearly the photoproduction limit, and Regge phenomenology; this 
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results in effectively a higher twist (or a 1/Q 2 type) contribution, that can reproduce the 
observed large, positive slope of the structure function ratio at small values of x, x < 0.05 
Q. As has been pointed out [|], a perturbative approach to the problem, addressed, similar 
to the free nucleon case, via the partonic content of the bound nucleon in the leading twist 
(LT) approximation, would result in small or essentially no Q 2 dependence of structure 
function ratios; such a partonic approach should also be more applicable at larger Q 2 . A 
perturbative approach, in principle, also includes a higher twist (HT) contribution, that dies 
away as 1/Q 2 - Such a contribution would therefore not be significant for Q 2 much larger 
than, say, 10 GeV 2 . There is no compulsive evidence for such HT contributions even down 
to 1 GeV 2 in the free nucleon case |7j. However, it is not known, at present, whether the 
Q 2 dependence of free and bound nucleon structure functions are the same. The current 
Sn/C data at small x, which show a significant Q 2 dependence, are precisely in the region 
of 1 < Q 2 < 10 GeV 2 where it is not clearly established that HT effects are subleading. It is 
pertinent, therefore, to ask what is the smallest Q 2 down to which a pure LT contribution 
to bound nucleon structure functions is still consistent with the new Sn/C measurement. 
The answer, surprisingly enough, is that a pure DGLAP approach is still consistent with all 
available data, within the experimental errors. We now present the details of the calculation, 
and substantiate our claim. 



2 Brief description of the model 

Specifically, we choose, for the purposes of effecting this comparison, a model |§ that uses 
standard leading order DGLAP || (or leading twist) evolution equations to obtain the 
nuclear structure functions at various Q 2 values, starting from a low input scale, Q 2 = 
fi 2 = 0.23 GeV 2 . The Gluck-Reya-Vogt (GRV) |?[] parametrisation is used for the free 
nucleon input, which is well-known to provide a good fit to the available proton data; the 
corresponding bound nucleon input densities are obtained by calculating distortions of the 
free ones due to the effects of nuclear binding and swelling. Binding, in particular, plays a 
special role at low-x, where shadowing is observed. Details of the model can be found in 
Ref. [|J. In summary, the increased radius of the bound nucleon, due to nucleon swelling 
Toll , results in a redistribution of the partons inside the nucleon ||11|| . The increase in the 



radius, 5a, of a bound nucleon, relative to a free one, is parametrised as, 

^ = Mi-i-i^ 1/3 ) ; (i) 

the free parameter, 5 vo \, was chosen to be 15%. Eq. (1), therefore, fixes the swelling for all 
A. Conservation of parton number and energy-momentum, along with use of Heisenberg's 
uncertainty principle, then specifies the input bound nucleon parton densities, starting from 
the free nucleon ones, using only one free parameter that parametrises the extent of swelling, 
and with fixed A dependence, as indicated in eq. (1). Nuclear binding is introduced by 
assigning the consequent energy loss, as computed using the Weizsacker mass formula, to 
the sea quark sector; the sea quark density is thus depleted in comparison to the free nucleon 
case, by an amount dependent on its mass number, A. Such a depletion also occurs at the 
time of scattering, for partons with small x, x < 0.1 f|. The model predicts structure 
functions with very few free parameters, and with well-defined A dependences. Results have 
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been obtained for various structure function ratios at average values of Q 2 , corresponding to 
the measured values; good agreement with data has been found in the region of validity of 
the model. (The model is not valid at large x, x^O.5, since Fermi motion effects have not 
been included). 



2.1 Comparison with data 

In view of the new Sn/C data, we now examine in detail the Q 2 dependence of the 118 Sn and 
12 C nuclear structure functions and their ratios within this dynamical model. We emphasise 
that, once the model has been fitted to data from nuclei such as Helium or Carbon, as has 
already been done, the resulting Tin structure function is essentially fixed, with no more 
free parameters. Furthermore, the magnitude and the Q 2 dependence in the base case, viz., 
for the structure function of proton or deuteron (which is considered to be the prototype 
"average nucleon") is well reproduced by this model. (It is in fact tailored to do so). Since 
the bound nucleon densities are computed in a manner similar to that of the free nucleon 
ones, that is, by the DGLAP evolution equations, the Q 2 dependence of the bound densities 
is similar to that of the free ones. The only difference lies in the input (starting) densities, 
which are computed according to standard nuclear inputs. Hence our model exhibits only a 
weak Q 2 dependence for all structure function ratios [I2~| . This can be seen in Fig. [I], where 



the model prediction for the structure function ratio, R = F^/F^, is plotted for various 
x values in the range between about 0.01 and 0.4 (the central value of the corresponding 
experimental x bin was used) el. S 9 function of Q 2 , for 1 < Q 2 < 100 GeV 2 . Note that all the 
data lie in the accepted perturbative regime of Q 2 > lGeV 2 . 

Although the central values of the data at low-x are strongly increasing with Q 2 , it is 
seen that the model is in good agreement with the available data, within the errors of the 
experiment. The \ 2 P er degree of freedom of the fits, as shown in Table [I], is seen to be 
around unity or less. The only exception is the last x-bin, x = 0.55; however, here our model 
is not expected to hold, since we have not included Fermi motion effects. (This also seems 
to indicate that these effects are severely A-dependent and so do not cancel in the ratio of 
structure functions). The average x 2 /d.o.f. = 138/142 for the complete data set, indicating 
that the model is consistent with the data. This further improves to 109/136 if we drop the 
last data point, x = 0.55. 



X 


X 2 /d.o.f. 


d.o.f. 


X 


X 2 /d.o.f. 


d.o.f. 


X 


X 2 /d.o.f. 


d.o.f. 


0.0125 


0.59 


8 


0.055 


0.78 


11 


0.25 


0.53 


11 


0.0175 


1.18 


9 


0.07 


0.93 


12 


0.35 


1.10 


9 


0.025 


0.65 


9 


0.09 


0.27 


12 


0.45 


1.16 


7 


0.035 


0.96 


11 


0.125 


0.90 


14 


0.55 


4.77 


6 


0.045 


0.55 


11 


0.175 


0.99 


12 









Table 1: The \ 2 fits °f the model with the NMC data on the Q 2 dependence of the Tin to Carbon 
structure function ratio [|| , for various x bins, the central values of which are shown. 



Hence, we see that a purely DGLAP-type perturbative evolution of an input set of densi- 
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ties predicts a very small Q 2 dependence of the Sn/C structure function ratio ||, |5| . Moreover, 
the prediction is consistent with the available NMC data. 

The Q 2 dependence of the structure function ratios, R A = F^/F®, for A = Li, C, Ca, 
and Pb, also show very little Q 2 dependence in the shadowing regime, x < 0.1. The resulting 
fit to a linear Q 2 dependence, of the form, 

^ = Mn Q 2 + c, (2) 

yields a slope, b, that, though positive at small-x, and negative at large x (x > 0.1), is almost 
vanishing throughout, and compatible with the corresponding measurements available from 
NMC and E665 |l], |2|] in the region of overlap (we do not include the points which have 
x > 0.5 and Q 2 < 1 GeV 2 since we are using a perturbative evolution equation without 
Fermi motion effects). We therefore believe that the model predictions are compatible with 
the available data. 

This implies that the Q 2 dependence of the bound and free structure functions are essen- 
tially the same. Note that the corresponding free nucleon data in the same x region can be 
fitted with a DGLAP type Q 2 behaviour and do not seem to need inclusion of HT effects in 
the region 1 < Q 2 < 10 GeV 2 . Hence the Q 2 dependence of free and bound nucleon structure 
functions seems to show the same origin, viz., parton splitting. We test this by comparing 
the Q 2 dependence of the Carbon structure function data alone, not just structure function 
ratios, against the preliminary high statistics data available from the NMC |13|]. This is 
shown in Fig. 0, for the same ranges of x, and for the same Q 2 bins as the Sn/C data. Note 
that the error bars shown here correspond to statistical errors alone. We see that there is 
reasonable agreement with the data for all x ranges. However, the x 2 /d.o.f. here is worse, 
being only 122/87; this is because of worse fits for 0.04 < x < 0.1 as can be seen from 
Table |^. This can be greatly improved by an overall normalisation change and does not 
require a change in the slopeQ For instance, a decrease in the overall normalisation by just 
3% decreases the x 2 /d.o.f. to 74/87 (or to 68/83, if we drop the last data point), which is 
very acceptable. However, we do not do this here, and merely keep in mind that there would 
also be an improved agreement on using a next-to-leading (NLO) order fit rather than just 
a leading order one. 



2.2 The polarised case 

Finally, we consider the Q 2 dependence of the corresponding polarised structure function, 
g^(x,Q 2 ). In particular, the structure function of interest |TJ| is that of 3 He, since it was 
the target used in determining the neutron spin dependent structure function |T5| . We use 



an extension of the same dynamical model to the polarised case [pi 16 1 to study the Q 
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dependence, where the free nucleon input now corresponds to the fits from the Gluck-Reya- 



lr Tliis can be accomplished by a small change in the extent of swelling, which remains a free parameter in 
this model. This would cause a corresponding change in the structure functions of all other nuclei; however, 
this should not worsen the agreement of the model with these other measured structure functions. 
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Table 2: The x 2 fits of the model with the NMC data on the Q 2 dependence of the Carbon structure 
function, F 2 c (x,Q 2 ) pi, for various available x bins, the central values of which are shown. 



Vogelsang (GRVs) [JIT] parametrisation. We see from Fig. g that the polarised ratio, 

He/n 

£=^V, (3) 
9i 

is more sensitive to Q 2 than the corresponding unpolarised ratio, 

pHe/n 

R n = -^r , (4) 

where the superscript, n, on the unpolarised ratio indicates that the "average neutron" 
structure function in He was computed, in order to effect a meaningful comparison with R. 
(See Table |3| for a comparison with data.) At small x, where the sea densities dominate, 
the Q 2 behaviour of R and R are different, owing to the different Q 2 dependence of the 
corresponding unpolarised and polarised splitting functions, Py and AP^-. At larger x, from 
x>0.1, since P^ q = AP qq , i.e., the nonsinglet unpolarised and polarised splitting functions 
are the same, R and R show the same Q 2 behaviour. 

Such a Q 2 dependence [|I7]] implies that the corresponding spin asymmetry, 



A/n 

which is what is usually measured in DIS with polarised targets, will show more sensitivity 
to Q 2 . This can be seen in Fig. |4] where both the free (A™) and bound (A^ e ^ n ) nucleon 
asymmetries are shown (using F 2 = 2xFi) for various x values accessible by current ex- 
periments. The extent of the Q 2 dependence, as has been exhibited by a comparison with 
the available data in Table [3], gives hope that, when more data is available, this large Q 2 
dependence can actually be observed. For instance, the asymmetry changes by about 15- 
20% for x from 0.025-0.2, in going from Q 2 = 1 to Q 2 = 10 GeV 2 . As has already been 



pointed out |16|], the nuclear dependence effectively cancels in the ratio, eq. (5), so that the 
free and bound nucleon asymmetries are nearly equal; hence the E142 measurement (which 
measures the bound nucleon neutron asymmetry in Helium) still provides a good measure 
of the free neutron asymmetry, and hence no information at all about nuclear effects. It 
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Table 3: The model prediction for the polarisation asymmetry in Helium compared with data 
[15] at the same central values of (x,Q 2 ) as the data; errors shown are statistical and systematic, 
combined. The Q 2 dependence of the computed asymmetry is indicated by showing the asymmetry 
at Q 2 = 1, 10 GeV 2 in the last column. 



seems as if information on spin dependent nuclear structure functions can only be obtained 
by measuring ratios similar to that shown in Eq. (3), or plotted in Fig. [| 

rather than from measurement of the conventional polarisation asymmetry, A^ n . 



3 Discussion and Conclusion 

We conclude that a perturbative approach to nuclear structure functions with the use of 
the leading twist DGLAP evolution equations predicts, not just the ratios, but also the 
magnitudes of various nuclear structure functions, and their Q 2 behaviour as well, that are 
consistent with all data so far available. Such an approach typically yields a very flat Q 2 
behaviour for structure function ratios at small x. A similar Q 2 behaviour of bound and free 
nucleon spin independent structure functions is therefore indicated. Another perturbative 
approach, using parton recombination and Q 2 rescaling arguments, has been able to obtain a 
non-trivial Q 2 dependence , similar in nature to that obtained by non-partonic approaches 
0, 0| . The Q 2 dependence in this model, though, is very sensitive to the input parameters 
0. However, a persistence of the significant Q 2 dependence in the Sn/C ratio at small x (for 
Q 2 > 1 GeV 2 ), accompanied by a shrinking of the corresponding error bars by the acquisition 
of more data, could well signify the presence of higher twist effects, which have not so far 
been observed for the free nucleon case at such Q 2 values. This would significantly affect 
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the extraction of spin dependent structure functions from polarised deep inelastic scattering 
experiments off polarised nuclear targets as well. Hence, a detailed experimental study of 
the Q 2 dependence of the spin polarisation asymmetry is also urged. In short, it would be 
interesting to obtain more data for various nuclei in the small-x region, both in polarised and 
unpolarised experiments, in order to verify if a universal Q 2 behaviour of free and bound 
nucleon structure functions exists. This is bound to throw more light on the nature and 
origin of the modifications of parton densities in nuclei. 
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Figure 1: The model prediction for the Q 2 dependence of the structure function ratio for Sn/C, in 
comparison with data from the NMC |J, with statistical and systematic errors added in quadrature. 
Average (central bin) values of x are shown. 
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Figure 2: The model prediction for the Q 2 dependence of the bound nucleon structure function in 



Carbon for various values of x. Data shown correspond to preliminary results from the NMC [13] 
with only statistical errors shown. 
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Figure 3: The model prediction for the Q 2 dependence of the spin independent (R) and spin 
dependent (R) structure function ratios, for He/n, are shown as solid and dashed lines respectively 
for various x values. 




Figure 4: The model prediction for the Q 2 dependence of the polarisation asymmetry, A^ e (dashed 
lines) and the free neutron asymmetry, A\ (solid lines), for various values of x from 0.0055 to 0.5. 
A comparison with available data |l5| can be found in Table [|. 
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